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Mass dependence of the hairpin vertex in quenched QCD 

Stephen R. Sharp^ 

Physics Department, University of Washington, Seattle, WA 98195-1560 

The pseudoscalar "hairpin" vertex (i.e. quark-disconnected vertex) plays a key role in quenched 
chiral perturbation theory. Direct calculations using lattice simulations find that it has a significant 
dependence on quark mass. I show that this mass dependence can be used to determine the quenched 
Gasser-Leutwyler constant L5. This complements the calculation of L5 using the mass dependence 
of the axial decay constant of the pion. In an appendix, I discuss power counting for quenched chiral 
perturbation theory and describe the particular scheme used in this paper. 



I. INTRODUCTION 



O 

00 

o 



X 



cn ■ 

o : 
o. 

' Although the quenched approximation to QCD is being gradually superseded by unquenched simulations, it remains 
^ ' useful as a qualitative guide to the physics of QCD, and as a testing ground for new methods of discrctizing light 
1^ , fermions. In the latter role, the unphysical chiral singularities of quenched QCD turn out to be helpful — reproducing 
^***| ■ them provides evidence that the light fermion method under study can approach the chiral limit. Examples of these 
singularities are the enhanced chiral logarithms in the pion mass and the pseudoscalar decay constant P, Q . The 
source of these singularities is the hairpin vertex, i.e. the appearance of a double-pole contribution in the flavor-singlet 
pion two-point function. In infinite volume (where the contributions of exact zero modes vanish) this vertex is the 

■ key unphysical feature of quenched QCD. 

J> \ Quenched chiral perturbation theory (QxPT) allows one to systematically predict the form of these singularities. 
0^ ' Although QxPT does not have the same theoretical underpinnings as chiral perturbation theory for unquenched QCD, 
and thus must be treated as a phenomenological theory, it has been quite successful in describing theproperties of 
light pseudo-Goldstone bosons. A particularly thorough study has been undertaken by Bardeen et al. 0)0] (referred 
to as BDET in the following).^ They find that consistent values of the hairpin vertex can be obtained both from a 
, direct calculation of the hairpin correlator and from fitting to the chiral logarithms predicted from loop contributions 
' to various quantities. 

"Y^ I There is, however, a puzzling aspect to the results for the hairpin vertex. On the one hand, the vertex is found to 
. have substantial dependence on quark mass. As I argued in Ref. |lll |. this can be explained by a non-zero value for 
I ' the momentum-dependent part of the hairpin (a$ ^ 0.5 in the notation explained below). On the other hand, BDET 
Oh, note that, at tree-level, one can also determine a$ from the coefficient of the sub-leading single-pole contribution to 
^ ' the hairpin correlator. They find that a$ determined in this way is much smaller than that determined from the mass 
dependence, and, in fact, is consistent with zero. 

The purpose of this note is to resolve the discrepancy between these two methods for determining a$. It turns 
out that, at next-to- leading-order (NLO) in quenched chiral perturbation theory, the hairpin vertex picks up a mass 
dependence proportional to L5 (the quenched analog of corresponding Gasser-Leutwyler constant), in addition to that 

■ proportional to a$. The coefficient of the single pole, however, is not affected at this order, and remains proportional 
to a$ alone. This allows a consistent description of the hairpin results: one should use the BDET method to determine 
a$, with a result close to zero, while the mass dependence of the hairpin vertex determines L5. This provides a check 
of determinations of L5 based on the mass dependence of the axial decay constant. 



II. DETAILS OF CALCULATION 



I consider quenched QCD with Ny valence quarks. For the following considerations I need Ny > 2, and for 
definiteness I use Ny = 3. I take all the quarks to be degenerate, since this is what has been done in numerical 
evaluations of the hairpin vertex. The leading order quenched chiral Lagrangian is then Q,Q 

Co = ^Str {d^Ud^U^ -xU- U^x) + "^o*o + «$(5m$o)' , (1) 



'Electronic address: sharpe@phys. Washin gton . edu| 
1 See also the work of Refs. |lBI3IEIlEl- 
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where "Str" indicates supertrace, U = exp(2j<I>//) is an element of C/(3|3), $ contains the pseudo-Goldstone bosons 
and fermions, and $o = Str($)/-\/3 is the "super-ry"' field. Quark masses are contained within x = ^BqM, where 
M — diag(m, m, m, m, to, to) is the mass matrix of quarks and ghost-quarks. I assume that to is small enough that it 
is appropriate to use chiral perturbation theory to describe the properties of the pseudo-Goldstone hadrons. 

The parameters / (pion axial decay constant in the chiral limit, normalized so it is close to 90 MeV) and Bq are the 
quenched analogs of the corresponding parameters that appear in the unquenched chiral Lagrangian. The parameters 
Too and a$ are, however, special to quenched chiral perturbation theory, in which one cannot integrate out the <I>o. 
In the large Nc and chiral limits, one can show that toq = to,,' 1, 2]. Assuming this holds also at Nc = 3, one expects 
Too ~ 0.85 GeV. The presence of the new mass-scale toq gives rise to the enhanced singularities of the quenched 
theory, the strength of which is parameterized by 6 = TOq/ (487r^/^). It will be useful in the following to keep in mind 
the approximate value of S. The naive estimate (using QCD values for the parameters) is (5 ~ 0.2, while simulations 
find values in the range S = 0.1 — 0.2 (for a recent review see Ref. 22)- What is most important about these values 
is that they are small. 

The present calculation will be performed at NLO in a power counting scheme for quenched chiral perturbation 
theory described in the appendix. It will require the standard NLO Lagrangian of the unquenched theory, suitably 
generalized by changing traces to supertraces, plus an additional term: 

CiA - -Li [Str (a,,{/a^C/t)] ' - L2 Str {d^,Ud,U^) Str {d,,Ud,U^) - L3 Str {df^Ud^^U^ d,Ud,U^) 
+ L4 Str (9L/aC/t) Str {xU^ + Ux) + ^5 Str {dUdU^ [xU^ + Ux])- Lq [Str {xV^ + Ux)] ^ 

- Lj [Str (xC/t - Ux)]^ - Ls {Str (UxUx) + Str (xU^xU^) } 

- Lq Stv {d^U d,U^ d^U d,U^) . (2) 

Here Li_8 are the quenched analogs of the corresponding Gasser-Leutwyler coefficients. Lq multiplies an operator 
which is not independent in unquenched SU{2) or SU{S) chiral perturbation theory, but which is independent in 
the quenched and partially quenched theories. Its presence in partially quenched theories was noted and explained 
recently |13|. That it is also needed in the quenched theory can be seen by generalizing the arguments given there. 

In the quenched theory, one must also include other operators obtained from those in Cq and £1^1 by multiplying 
by powers of ^o- These introduce many new constants. It turns out, however, that the resulting operators do not 
contribute at NLO to the quantities I consider here, as explained in the appendix. For illustration, I list two of the 
leading operators of this type: 

CiB = «o(TO2//2)$4_^gig(,$^/j){Str(t/xC/x)-Str(xC/W)} , (3) 
where wo,8 are additional couplings. 

Following BDET, I consider the following correlation functions: 

Ac(p) = ^e-'f-^(QA75<ZB(x) ^759^(0)), (4) 
A„(p) - / e-'P--{qAl5qA{x) qBl5qB{0)) . (5) 

J X 

Here qA and qB are two different valence quarks, and the expectation value indicates the quenched average. The 
subscripts c and h are a reminder that only quark-connected contractions contribute to Ac, whereas A/i contains only 
quark-disconnected contractions. In other words, Ac is the usual non-singlet pseudoscalar two-point correlator, while 
A/j is the pseudoscalar hairpin correlator. These are the correlators used in simulations to study the hairpin vertex. 
At low momenta the connected correlator is well described by a single pole 

Ac (p) = - 2 + non-pole , (6) 

where /p is, by definition, the pseudoscalar decay constant, and the sign arises from Wick contractions. By contrast, 
the hairpin correlator is observed to have a leading double-pole contribution: 

where /p is to be determined from eq. ©. I will refer to the residues T) and S as the hairpin and single-pole 
coefficients, respectively. They are defined to be independent of p^, but can depend on to^, as do fp and M^r. 
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With this notation in hand, I can now state more precisely the observations made in BDET and previous work. 
These are, first, that 5 « 0, i.e. that A^(p) is consistent with a pure double-pole; and, second, that V has a significant 
mass dependence. I now turn to the predictions of the quenched chiral Lagrangian for 2? and S. 

To begin with I need to determine the operator in the effective theory which matches onto the pseudoscalar density. 
This can be done, as in unquenched chiral perturbation theory, by taking appropriate derivatives with respect to the 
quark mass matrix, treating the latter as a source. At leading order, using £o, this gives qA^sQB 2iBaf^AB+0{^^), 
for all A and B, leading to the tree- level results 

Acb) = -^¥r^ (8) 

P +x 

^"-^^ ^ ^ (7T^- 

Comparing to eqs. and {Tj) I find 

Ml=x; fP = 2B,f; 5 = f . (10) 

These leading order results are well known and agree with those quoted in BDET when normalization and sign 
conventions are taken into account. 

As noted above, at this order a$ provides both the mass dependence of V and the value for S. It is perhaps worth 
understanding this point in more detail. One might have thought that a term in Co of the form 

{f/A)v2Str{xU^ - Ux)t^o/f (11) 

would provide additional mass dependence to the hairpin vertex, and thus break the correlation between V and S. 
However, this term also contributes to the effective field theory expression for the pseudoscalar density, and leads to an 
additional term in S, in such a way that V2 can be absorbed into a redefinition of a$. This is a consequence of the fact, 
noted in Ref. 0, that the V2 term ((TT|l can be removed by a field redefinition [expficitly this isU ^ U ex.p{—iv2^o/f)] 
and a redefinition of parameters (a$ — > — V3v2). 

Before proceeding to the details of the NLO calculation, it is necessary to discuss the power counting appropriate 
to quenched chiral perturbation theory. I give an overview here and a comprehensive discussion in the appendix. 
Standard unquenched chiral perturbation theory is, in the meson sector, an expansion in ~ M^/A^ ~ p^/A^, 
where A^ — 47r/. The issue in quenched chiral perturbation theory is how to treat the extra constants that appear, 
the most important of which is toq- This appears in two dimensionless combinations in the calculations considered 
here. First, loop contributions are suppressed by powers of S = TOo/(3A^), which, as noted above, is numerically 
small. ^ Second, the ratio of the two terms contributing to V in eq. (|10|l is a,j>M^/mQ. It is also reasonable to treat 
this combination as small, given that nio is close to A-^ and a$ is small. To do so consistently, however, one must 
develop a systematic power counting scheme in which both combinations are small, and for which one can assess the 
relative size of other contributions. 

Two such schemes are explained in the appendix. In the first (the standard or PCI scheme) one takes 5 ~ ag, /3 ~ e^, 
while in the second (the PC2 scheme) S ^ ^ e^^"^. In both schemes the q;$ contribution to V is non- leading: 



2 



rUn 3 S 



,2 



e 



PC1/PC2. (12) 



The schemes differ, however, in the size of other contributions. In the PC2 scheme, since 6 is larger than e^, "second- 
order" terms of size S'^, Se^, and are progressively smaller, while all three are treated as of the same size in PCI. 
Thus the PC2 scheme has some similarities to that used in BDET, in which terms of relative size Se'^ log(M^//i) are 
kept, but those of size e'^ are not. This point is discussed further in the appendix. 



^ In fact, the leading logarithmic corrections are proportional to powers of 5 log(Af7r//i), with fi the regularization scale, and thus dominate 
over the leading term in the extreme chiral limit. For most practical simulations, however, including those of BDET, the logarithm is 
of 0(1), and does not lead to a significant enhancement. I assume this to be case in the power counting schemes I use. 
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FIG. 1: Examples of diagrams contributing to the pseudoscalar hairpin correlator which are associated with vertex renormal- 
ization. Simple vertices are from the first two terms in Co, crosses indicate mo or a# vertices, filled boxes are interactions 
from C\A, with the specific term which contributes being noted. Note that if the rightmost propagator and hairpin vertex are 
removed, these become contributions to the quark-connected correlator. 

In this paper I calculate the connected and disconnected correlators at NLO in the PC2 scheme. By NLO I mean 
that corrections of size up to and including relative to the leading contribution are included. After amputation, the 
disconnected vertex has the form V + {p^ + M^^)S + . . ., from which one can extract V and S. Note that one obtains 
a lower order result for S than for V because of the explicit factor of {p^ + M^) multiplying S. 

I choose the PC2 scheme because it reduces the number of diagrams that contribute and thus simplifies the 
calculation. In both schemes, however, the general point is the same: there are other terms which contribute at the 
same order as the a$ vertex, and these must be included. 

The calculation of the connected correlator is standard and leads to the replacement of fp and in eq. lO with 
their NLO expressions. The leading diagrams that contribute can be deduced from those in figs.nandiniby removing 
the rightmost propagator and the hairpin vertex to which it couples. As can be seen from tabled a NLO calculation 
involves only diagrams (a) and (b) from both figures — other diagrams are of higher order than in both power 
counting schemes. The explicit expressions for fp and are given in BDET and elsewhere, with the latter repeated 
in the Appendix, eq. ljA8|l . They are not needed to determine I? and S. 

For the disconnected correlator, one can divide the contributions into the four classes. (1) Those renormalizing 
the external pseudoscalar density, and which would be present also if the density were flavor non-singlet. Examples 
of these are shown in fig. ^ (2) Those renormalizing the external density, but which are present only because the 
density is flavor singlet. These contribute only to the single pole amplitude S, but not to V. Examples are shown in 
fig- 13 (3) Those renormalizing the mass and wavefunction of the pseudo-Goldstone boson, shown in fig. 01 (4) Those 
associated with the hairpin vertex, shown in fig. ^ Two observations greatly simplify the calculation. First, the 
contributions of classes (1) and (3) are identical to those appearing in the connected correlator. This means that the 
same expressions for fp and apply to both connected and disconnected correlators — as one would have expected. 
Second, all contributions of classes (2) and (4), aside from the leading order diagram fig.EJa), are smaller than NLO."^ 

These observations imply that the sole NLO correction to the hairpin vertex in the PC2 scheme comes from 
wavefunction renormalization. This in turn arises at NLO only from fig. |2Ib), as fig. Ola) only renormalizes the 
mass. Since there are two boson propagators in the disconnected correlator, compared to only one in the connected 
correlator, only "half" of the wavefunction renormalization in the disconnected case is absorbed into the one-loop 



^ This simplification does not occur in the PCI scheme, in which at NLO one must also calculate figs. I^a-b*) andHTb-cV 
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FIG. 2: Examples of diagrams contributing to the pseudoscalar hairpin correlator which contribute directly to the single pole 
residue S. Notation is as in figs. except that boxes with attached arrowheads are from Cm (the arrowhead indicating the 
$0 field). 




FIG. 3: Examples of diagrams contributing to the pseudoscalar hairpin correlator which are associated with pion mass and 
wavefunction renormalization. Notation is as in fig.Q Note that if the rightmost propagator and hairpin vertex are removed, 
these are contributions to the quark-connected correlator. 

result for fp. The other half renormalizes the hairpin vertex, leading to the results^ 

2? = i|m^(l-8L5-^)-aa4^}; S = ^. (13) 



* In an earlier version of this paper I also kept corrections of relative size eS^ log(M.n-//i), coming from fig. Ole), as in the calculation of 
BDET. The analysis in the appendix shows that in either power counting scheme one should then also include all the other diagrams in 
figs. |3]and|l] including non-trivial two loop diagrams. I have not attempted this calculation. 
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FIG. 4: Examples of diagrams contributing to the pseudoscalar hairpin correlator which are associated with renormalization 
of the hairpin vertex itself. Notation is as in figs. Q and |21 

Note that the NLO contributions do not change S from its tree level value in eq. 11U|) . In particular, the correction to 
S from wavefunction renormalization is a correction of relative size e"* in the amputated disconnected correlator, and 
thus beyond the order I am considering. Note also that one can interchange x and in the result, as the difference 
is of higher order than I consider. 

Equation is the new result of this paper. It shows that the mass dependence of the hairpin T> is not tied 
to the coefficient of the single pole S. Thus, in particular, the single-pole term can be absent (as found, to good 
approximation, by BDET), but the hairpin vertex can have a linear dependence on M^. Furthermore, the latter 
dependence provides an alternative method for determining L^. 

There is a nice consistency check on the result (I13II . At p = 0, the hairpin correlator is proportional to the topological 
susceptibility of the pure gauge theory 

m^Ahip^O)^Xt^{iy^)/V, (14) 

where m is the quark mass, v the topological charge, and V the four dimensional volume. This is the "fermionic" 
method for measuring the susceptibility. Since this formula holds for any m, and Xt, being a pure gauge quantity, is 
independent of m, the left-hand-side should also be independent of m. Inserting cq. H13() into the general form lO I 
find 

A.(P = 0) = ^ + SM'^ + 0{Mt)) (15) 



3M4 V P 
= 4?(l-8L5-^) . (17) 

where in the last step I have used the axial Ward identity 2m fp — M^/a, which holds to all orders in chiral 
perturbation theory (and can be checked at NLO from the expressions in BDET).'^ Using the expression given in 
BDET for the axial decay constant /a, I find 

A{p = 0) = ^. (18) 



^ Here I use the BDET normalization such that fj^ = V^f in the chiral limit. 
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This is independent of the quark mass, as required. 

One can try turning this argument around, i.e. assuming H18|l and deriving eq. H13|) . This will, however, only 
constrain the combination T) + M^S, and not allow a separate determination of V and S. 

I close this section with a technical comment. In order to predict the coefficient of the single pole in the hairpin 
correlator one must control the corrections to the operator which creates the meson at 0{p^) relative to the leading 
order. This is because such a correction multiplying the double-pole leads to a contribution to the single-pole. Thus 
it is crucial that one calculate the hairpin correlator with an operator which can be matched onto the effective theory 
without introducing additional coefficients at 0{p^). This is only true for the pseudoscalar density, used here, and 
for the axial current. With any other operator coupling to the pseudoscalar meson, e.g. qd^j^q, the coefficient of the 
single pole could not be predicted. 

In light of this discussion, one might be concerned about the need to include the operator 

Str(x52[C/ + [/t]) (19) 

in CiA, since this would give a contribution to the pseudoscalar density proportional to p^$. It turns out, however, 
that this operator can be removed (as in unquenched chiral perturbation theory) by a field definition and appropriate 
shifts in L5, Lg and a$. Thus it is a redundant operator.^ I have checked this by explicitly calculating its contributions 
to /a, fp and Ac,h- 

III. CONCLUSION 

In this note I have presented a complete NLO calculation of the hairpin vertex in quenched chiral perturbation 
theory. The result, eg. Ijl^ll , resolves a puzzle concerning the determination of the parameter a$ . The method I had 
proposed in Ref. [Tl|, using the quark mass dependence of the residue of the double-pole, turns out to be incorrect. 
The correct method is to use the residue of the single-pole, which, according to BDET, leads to « a$ ^ 1. The 
mass dependence of the double-pole residue then allows one to determine L5 . 

It will be interesting to reanalyze the results of BDET and earlier work in this light, and compare the determinations 
of L5 with those obtained from the mass dependence of the axial decay constant, /a- As a rough first attempt, I note 
that my earlier estimate a$ « 0.5 converts to ^5 « 0.5/(487r^^) « 1.3 x 10^'^, where I have taken 6 = 0.1 from BDET 
for the value of the hairpin vertex in the chiral limit. This is close to the values for L5 quoted in BDET. 

The calculation also suggests a simple way of displaying the consistency of the different determinations of the 
parameters of quenched chiral perturbation theory. It follows from the previous considerations that 

{V + M'^S)fl = Axt + OiM^). (20) 

In the notation of BDET, this means that the combination [(ttiq**)^ + M^a^] x should be independent of m at 
linear order. Since a$ is small, this constraint explains the observed result that mg^ decreases with increasing m in 
terms of the fact that /a increases with m. 

When applying the results of this paper to numerical simulations, one should keep in mind that the previous 
discussion did not account for finite volume or discretization errors. Concerning the former, one can include the 
leading effect by replacing the loop integrals with those appropriate to finite volume. As for the discretization errors, 
if one uses Wilson-like fermions, then there will be 0{am) corrections to the results presented here, unless the action 
and operators are non-perturbatively 0{a) improved (which is not the case in the work of BDET). This is because of 
the breaking of chiral symmetry by the action. In particular, the quantity in eq. I|20|) will have a mass dependence 
proportional to am. Thus a possible use of this relation is to gauge the size of such discretization errors. 
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APPENDIX A: POWER COUNTING IN QUENCHED CHIRAL PERTURBATION THEORY 

In this appendix I discuss power counting schemes for quenched chiral perturbation theory, and explain that used 
for the calculation described in the main text. I also take the opportunity to collect some results that are presumably 
known to workers in the field but which, as far as I know, are not available in the literature. I use the standard 
notation for the chiral power counting parameter: — x/A'^, where x refers to either or p^. 

To make the arguments I will frequently use the fact that, in the quenched theory, diagrams involving pions (and 
their graded partners) can be converted in a precise, well-defined way, into those involving quark (and ghost quark) 
lines. This involves followingflavor indices, and is allowed because of the necessary inclusion of the $o field in quenched 
chiral perturbation theory ]^ Il7| . The advantage of the quark-line representation is that the impact of quenching 
is transparent: all diagrams with internal quark lines cancel against those with corresponding ghost quark lines. In 
the following I only consider the diagrams which remain after such cancellations have occurred, i.e. I consider power 
counting only for diagrams that give non-trivial contributions. 



1. Power counting with £o 

The largest contributions to correlators involve vertices from Co, eq. Q), so I consider first the power counting 
appropriate for this Lagrangian alone. For diagrams with no TOq and a$ vertices, standard chiral power counting 
gives the following result for the amputated N-point correlation function with external momenta of 0(M^):^ 

An ^ fN-2 (^^)"'^ 1 ""^t — number of loops . (Al) 

For example, at leading order the two-point correlator gives the pion mass, Ai ^ % ~ M^, with corrections suppressed 
by the number of loops. 

This result, and those that follow, also apply to the N-point correlator of the pseudoscalar density, except that 
eq. IjAip must be multiplied by Since what matters is the relative size of different terms, and not the overall 
factor, the ensuing analysis will be directly applicable to the correlators considered in the main text. 

Including TOq and vertices, I note that each such vertex introduces an extra propagator compared to the diagram 
without the vertex, and so their contributions (after subtracting divergences so that loop momenta are of order M^) 
are 

mo vertex ^ ~ ^ = ^ ' vertex ^ - a*/3 . (A2) 

I include the factors of 1/3 since, given the definition $o = Str((f>)/-\/3, they are generic. Including these results the 
power counting for Co is 

AN-j^ieT-'''(Sr{^y\ (A3) 

where ns and Ua are, respectively, the number of itiq and a$ vertices. The appearance of in the contribution 
of the mo vertex in eq. ljA2p is the reason for the modification of power counting in quenched chiral perturbation 
theory. As discussed below, one typically chooses 6 > e^, so that additional itiq vertices either come at "no cost" or 
are enhanced. 

There is, however, a constraint on the number of TOq and a$ vertices due to the cancellation of internal quark and 
ghost loops. A simple example of this constraint is that one cannot have two such vertices on the same propagator, 
since this would require an internal quark loop. As I will explain, the general form of this constraint is 

ns + na <ni + ricut , (A4) 

where ricut is the number of independent purely gluonic "cuts" through the correlator. In other words, if one traces 
the quark lines, they will fall into ncut + ^ quark-disconnected components. In perturbation theory the only connection 
between these components is provided by gluons. In the cases of interest here, the connected correlator Ac{p) has 
TT-cut — 0, while the hairpin correlator Ahip) has ricut = 1- 



This also holds for the corresponding IPI correlator, and thus for the N-point scattering amplitude, although I will not use this here. 
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To demonstrate the result (jA4|l I first consider correlators with ricut = 0, i.e. those which are connected in terms 
of quark lines. I need the result that the four and higher-point vertices coming from Cq involve a single supertrace, 
and thus are themselves singly connected in terms of quark line diagrams. The cancellation of diagrams with internal 
quark loops then implies the following: there must be at least one path through the (pion language) diagram from 
each ttiq or a$ vertex to each of the external pions which does not pass through another rriQ or a$ vertex. This 
must be true if one leaves the original vertex in either direction. Were this not the case, there would be a segment of 
the diagram completely "cut-off" by TOq and q;$ vertices, corresponding to an internal quark loop in the quark-line 
language. 

A corollary of this result is that, if one were to "pull apart" each of the m§ and vertices to give two extra 
external pions for each such vertex, then the resulting "expanded" pion-language diagram is singly connected (in the 
usual sense). 

Now imagine starting with an unquenched chiral perturbation theory diagram and adding as many TOq and q;$ 
vertices as possible on the propagators. When this has been done, the corresponding expanded diagram must be a 
tree diagram, since if there were any loops, the internal propagator(s) could be expanded by adding additional rriQ or 
vertices. It must also be a singly connected (pion-language) diagram from the corollary of the previous paragraph. 
Now, when one now sews the ttiq and a$ vertices back together again, each leads to an independent loop momentum 
in the original diagram, and thus ng + Ua = ni. Since this is the maximum number of such vertices that can be added, 
one finds the inequality ljA4|l with ricut — 0. 

The generalization to ncut > is straightforward. One can now have additional tjiq and/or a$ vertices, since ricut 
of them can be used to separate the ricut + 1 quark-disconnected components of the diagram. Thus one finds the 
general inequality of eq. ljA4|l . 

When using this formula two additional results should be borne in mind. First, there is, in general, no non-trivial 
lower limit on the number of ttiq and a$ vertices. One might have expected the constraint ns + Ua > ricut, i.e. that 
rriQ and a$ vertices are necessary to lead to quark-disconnected components. This is not the case, however, because 
loop diagrams involving connected vertices can lead to quark-disconnected correlators, as discussed in Refs. 0, 0| 
for the hairpin vertex. Second, even if a given choice of ns, and satisfies the constraint in eq. ljA4|l . there may 
not be a diagram contributing to An- Examples of this will be seen below. 



2. Including C\ and higher order Lagrangians 

It is straightforward to include the "standard" O(e^) Lagrangian, Cia, in the power counting. Each vertex from 
CiA gives an additional factor of e^, as in the unquenched theory, and the constraint on rig -f ria remains unchanged. 
It should be borne in mind, however, that since some of the vertices in £1^1 contain two supertraces, it becomes more 
complicated to determine the "quark-line connectivity" of a given graph. This has no impact on the generic power 
counting, but means that a larger number of diagrams will vanish in the quenched theory. 

As noted in the text, one must also consider operators which are built upon those in £0 a-nd Cia by multiplying by 
additional factors of $o//- Two examples are given in eq. (jSJ. The key observation is that these factors have the same 
power of 1// as vertices arising from the expansion oi U = exp(2i$//). Thus the power counting is unchanged, as 
long as ng and ria are generalized to be the number of vertices built upon the original ttiq and q;$ vertices, respectively. 
Thus, for example, ng becomes the total number of vertices of the form {^oY with g = 2, 4, 6, . . .. What is changed is 
the connectivity of the diagram. Each extra factor of ^0/ f requires an additional loop, or gluonic cut, in the diagram, 
since $0 closes off quark lines. Thus the loop constraint generalizes to ng + ria + < ni + ricut, where is the total 
number of extra factors of $0 in the diagram. Note that a ($0)^ vertex contributes both to ng (which it increments 
by one) and to (which it increments by two). 

Each vertex involving extra factors of $0 comes with a new coupling, e.g. z;o,8 in £ib. These couplings are 
suppressed in the large Nc limit because a $0 vertex is connected only by gluons to the rest of the diagram. Each 
additional $0// leads to a suppression by l/Nc- For example, if the quark loop is connected by two gluons to the 
rest of the diagram, there is an additional color loop, and the counting is g'^Nc = {g^Nc)^ /Nc ^ l/Nc- (The \/lVc 
contained in / is exactly that needed to create a state which is correctly normalized in the large Nc limit.) Thus I 
find that vg ~ l/Nc and vq ~ ^/N^- For a general diagram, the additional suppression is by {Nc)~'^'t' ■ 

I should stress that I am not claiming to extract the full Nc dependence of diagrams. There is also implicit Nc 
dependence in /, in some of the Gasser-Leutwyler coefficients, Li, and in TOq and a*. Indeed the latter two are both 
proportional to l/Nc 1], My idea is that, based on the general success of the OZI rule, it is likely that constants such 
as vg, and are suppressed relative to unity. Keeping track of the powers of l/Nc associated with these constants 
allows one to implement this as one deems appropriate. 

Higher order terms in the quenched chiral Lagrangian involve more derivatives (acting on either U or <I>o) and more 
factors of These can be accounted for, along with those from vertices in Lia, by introducing the quantity rii, which 
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counts the number of factors of x plus half the number of derivatives beyond those that would appear if the vertices 
were from the lowest order Lagrangian Cq. Thus a vertex from the standard chiral Lagrangian with six derivatives 
and one factor of x would increment ni by three. Overall, diagrams are multiplied by an extra factor of (e'^)"^. 

There is an ambiguity in this definition of ni: does, for example, the vertex proportional to (9<I>o)* count as ns — 1, 
= 2 and ni — 2 or as Ua — 1, — 2 and ni = 1? In other words, does it contribute a correction proportional 
to Se'^/N^ or to (a$/3)e"'/iV^? These are not the same size in some power counting schemes, e.g. those discussed 
below. The ambiguity arises because dimensions can be balanced either by powers of toq or powers of /. In fact, this 
ambiguity also arises for the vertices without derivatives, and I have chosen there the overall factor of Tig//'^"^ 
rather than 1/ f^^^. Fortunately, for the calculation considered in the main text, I can leave this ambiguity unresolved, 
as the operators in question do not contribute until higher order than that I consider. 

The preceding considerations together imply the following final power counting result: 

- jfe (,2)».-».+"i (^)n. ^^y- (^_^y' ns + n^ + n,<n, + n,„, . (A5) 



3. Applying the power counting 

I now apply the result (|A5|) to the examples considered in the text, namely the two-point quark-connected and 
disconnected correlators (i.e. A2 with Ucut — and 1 respectively). To do so I must decide how to treat 5, q;$ and 
1/Nc relative to each other and to e^. This depends, of course, on the size of quark masses considered. I have in mind 
mesons somewhat lighter than the kaon mass, so that « 0.2. 

I consider two possibilities for the relative size of the different terms. The first is 

PCI : ^ An^ (e ) , (A6) 

which one might call "standard" power counting. The key choices here are that S ^ e^, so that there is no extra 
enhancement of the TOq vertices, and that q;$/3 is small. One might argue that a$/3 ~ e, rather than e^, because it 
is proportional to l/Nc- I prefer the smaller choice given the numerical evidence that q;$ <C 1. 
The second scheme is 

r>r<0 ■ A " ^4/3 A X ^ 2\"g+"i-"a/3+2"o/3+n^/3 ,. , 

PC2 . d---^-e/ ^AT- (e ) . (A7) 

This enhances the contributions from $0 vertices relative to "standard" chiral vertices, and is the power counting 
used in the text. The power of e is chosen for a reason that will be explained below. In fact, most features of the 
power counting, including the application discussed in the main text, are unchanged if e^/'^ is replaced with e. 

Table HI collects the contributions to the quark-connected correlator which are of size or smaller relative to the 
leading order, in one or both power counting schemes. Table ^ gives the analogous contributions to the hairpin 
correlator, though only up to size e^°/'^. I list all contributions consistent with the constraint in eq. (jA5|) . but, as can 
be seen from the table, some are absent in quenched chiral perturbation theory. For example, the standard one-loop 
"tadpole' diagram (ni = 1, all other n's zero) requires an internal quark loop for ricut = 0, although, as observed in 
Refs. 10,01, it is present for Ucut = 1- Another constraint is that the only term with = 1 and ni = 0, namely 
the interaction quoted in eq. l(TT|l . can be (and I assume has been) removed by a field redefinition. 

To illustrate the difference between the power counting schemes, and the general nature of the chiral expansion 
in quenched chiral perturbation theory, I discuss the application to the pseudo-Goldstone pion mass. This can be 
obtained from the amputated quark-connected vertex. In the PCI scheme, NLO contributions to the quark-connected 
correlator are from the Gasser-Leutwyler constants (entering at tree level) and the ttiq vertex (at one loop) . The result 
for the pion mass is 



^ = 1 - S[log{x/fi^) - 1] + 8{2Ls - Ls)^ , (A8) 
X f 



in dimensional regularization and the renormalization scheme of Ref. |lfij . The a$ vertex first contributes at NNLO, 
and is but one of many terms, including three-loop diagrams with three ttiq vertices! 

The PC2 scheme is in part an attempt to systematize the power counting scheme used in BDET. This involves 
keeping corrections of size Se^ log(M^), while those of size are dropped. Here I am not treating log(Mn-//Lt) as large, 
so I must take S larger than e^, which leads to a scheme like PC2. The choice 5 ~ e^^^ is made (as opposed to 5 ~ e, 
say) so that the three loop contribution is smaller than the Se^ correction. Nevertheless TableUshows that keeping 
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TABLE L Leading contributions to the quark-connected two-point correlator in quenched chiral perturbation theory. PC gives 
the size of the diagram relative to A2 = X using the general power-counting formula, eq. (IA5II . while PCI and PC2 refer to 
the two specific schemes discussed in the text. All contributions up to relative size in one or both power-counting schemes 
are included. The references to fig. |3] refer to the parts of the diagrams to the left of the rightmost hairpin vertex (which are 
thus quark-connected). 
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TABLE IL Leading contributions to the hairpin (quark-disconnected) two-point correlator in quenched chiral perturbation 
theory. All contributions up to e^"/^ relative to the leading order term in one or both power-counting schemes are included. 
Notation as in Table ^ The references to figs. H and |1] refer to the entire diagram. 
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the Se^ correction, which is of size e^*^/^, requires keeping many other terms. Most have been calculated, but not all. 
Those not yet calculated are the non-trivial two- loop diagrams involving two mg vertices [Fig. Hd)], and a subset of 
the contributions of size ~ 5e^ from Fig. He) ■ The latter receives contributions not only from the L5 and Lg vertices 
(calculated in BDET) but also from the L3 and Lq vertices (not yet calculated). 

One can go to higher order in 5 if one assumes that log(M7r/^) is large, and keeps all terms of size [i51og(M^/^)]". 
This requires calculating only the tadpole loops, the simplest examples being figs. H^) ^^'^ (§)> these can be 
summed to all orders . This may be required when working at the smallest quark masses used in present quenched 
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simulations |Ml^.ll0|. 

The results in Table im are used in the calculation in the text, and are discussed further there. I only note here 
that an important simplification is that fig. Ofa) does not contribute to V and S. This would have been the leading 
correction to T> and the leading term in S. It does not contribute because the tadpole loops involving the TOq vertex 
do not give rise to wavefunction renormalization. This follows in turn because the vertex coupling two $0 fields to 
two non-singlet pions comes only from the mass term Str(xJ7) + h.c. and not from the derivative term Str((??79[/).* 
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This simplification does not hold for the four derivative single supertrace vertices with couplings L3 and Lq, which contribute through 
fig. |3e), as can be seen in the table. 



